TITLE 

TYPE 
ENTRIES 



NUMBER OF WORDS 
TEMPORARY STORAGE 

PRESET PARAMETERS 

LOCATION CONTENTS 
3 OOF OOaF 



UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 

LIBRARY ROUTIN E F6 -239 
By A. T. Nordsieck 

Integration of a system of ordinary differential equations 

with automatic control of integration interval. (SADOI Only) 

Closed - no program parameters - entered with link in Q. 

Normal entry is at left of OL. Special resetting entry: 

If this routine is entered at the right of 1L it will 

reset itself into the standard as-read-in condition, 

using the present values of the S6, S7 parameters and 

with the truncation accumulating locations reset. Use 

this entry only if a new sequence of integrations, is to 

be started (different initial conditions or different 

parameters) without re-reading-in this routine; and then 

only for the first entry of the new sequence, otherwise 

much time will be wasted and accuracy will be lost. 

Special interval controlling entry: if we enter at left 

of 5L, truncation information is preserved but new values 

of parameters S6, S7 take effect. This entry is used when 

£ Q is altered during a sequence of integrations. See S6 

parameter. 
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a to a + 3n-l for. continuing temporary storage; 0, 1, 

a + 3n to a + 6n-l for non- continuing temporary storage. 

Locations 3 to 8 must contain the following parameters 

during input and operation of this routine. 



SIGNIFICANCE 

N(a+i), i = 0, 1, ... n-1, are the variables x 
(independent variable), y^ y^, ... y . Initial 
values placed here. Results appear here. Auxiliary 
subroutine reads from here, x must be in first location 
a. 
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LOCATION CONTENTS SIGNIFICANCE 

h OOF 00(a+n)F n > 2. n is the number of equations including the 

"~™ """~"~ dx 

independent variable equation — = 1- N(a+n+i) are 

the scaled derivatives 2~ f . . Auxiliary must supply 

i 

these quantities for i = 1, 2, ... n-1. It need not 

-m „ _-m 



supply 2 f Q = 2 . a + 6j& < 1024. 



OOF OQmF 



1 < m < 37- m is an integer «uch that 2~ |f | are all 

guaranteed < 1 throughout integration range, m is 
provided so that the y and f» may be scaled independently 
for efficient use of the a + 'i and a + n + i registers. 
NOTE: To change this parameter after program has been 

read in it is necessary to overwrite location 5 with new 

-39 
value of m»2 , overwrite location (a+n) with new value 

of (2~ - 2~ ) and, if this routine begins at p, over- 
write location (p+127) with new value 50(p+100) 10(38-m). 

-JL 

6 OOF 00(J! -m)F m + 1 <£ Q < 38. 2 is the largest interval (increment 

in x) which the user wishes to permit. Furthermore, no 

matter what interval is actually used, the points 

x 5= (integer) 2 ° are guaranteed not to be skipped over 

provided the initial value of x is such a point. 

NOTE: To change this parameter after program has tyeen read 

in, change location 6 and use the appropriate special entry. 

7 OOF OOeF [3A e] > Jl Q - m+1. Also [3/4 e] < 50-m approximately. 

Here [3A e] means the integral part of 3 A e. The integer 
e specifies the accuracy desired by the user: Upon inte- 
grating over a range A x = 1, the y. 's will be in error by: 
not more than about 1 in the e bit for "well-behaved" 



equations (2 



df i 



, =1); not more than about 1 in the 
[3A e]th bit for extremely "ill -behaved" equations 
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LOCATION 



CONTENTS 



OOF OObF 



DURATION 



(2 



-m 



df 



dx 



SIGNIFICANCE 



■zQ 

"^2 ); and intermediate amounts for 



intermediate cases . For ranges of integration 

other than A x = 1, the error is proportional to 

the range. If the above rules put the error "beyond 

the 39 bit, it will nevertheless be about 1 in the 

39 bit. The computation time depends on e as 2 e ' , 

roughly . 

NOTE: To change this parameter after program has 

been read in, change location 7 and use the appropriate 

special entry. 

At location b, left side, is the entry to the auxiliary 
subroutine , which must be supplied by the user and which 
computes the scaled derivatives 2~ m f. (i = 1, 2, ... 
n-l) from the y and places them in locations a + n + i. 

T = 1.20 [lO.^n + k t] ms, where T is average time for 
one useful integration step; t is time to execute 
auxiliary once. The factor 1.20 is accounted for by 
discarded integrations incidental to keeping the interval 
correct. When modified into a frozen- interval program 
functioning like code Fl, the time is reduced to: 
T' = 9-7 n + k t ms. per step. 



To Freeze Interval: 



If this routine starts at location p, 



Overwrite: 

p+12: L5 107L LO 78L 
p+i+6: L5 H^L hO k?L 
p+67: kOF L7F 



With: 

22 (p+2l)F OOF 
26 (p+3^)F OOF 
26 (p+69)F OOF 
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(interval will be fixed at the value 2 , no tests will be made, no integrations 
discarded, e will be ignored.) 

To Reverse Integration Sense: After read in, if this routine starts at p, 



Overwrite: 

p+38: k2 kkL L5 F 

jp+120: kO a+4n+l LO a+5n+l 

a+H: 2- m -2- 59 



With: 

k2(v+kk)F LI ( )F 

kO a+4n+l iA a+5n+l 
2 -m 2 -39 



FF Stops During Read-in: 



Order With 
Sexad. Address 

FF Okk 
FF 0^5 
FF Ok6 
FF 0^7 
FF 048 
FF 049 
FF 04K 



From Location: 

7L left 

9L left 

12L left 

15L right 

16L left 

l6L right 

19L right 



Significance: 

a + 6n too large 

n too small 

m too large 

X. too large 

m too small 

£ too small 

e too small 



FF Stops During Operation: 



Order With 
Sexad. Address 

FF QhS 



FF O^N 



From Location: 
53L left 

62L left 



Significance: 

A test, so mild that any number 
whatever should pass it, is failing, 
hence program is out of order. 

The interval must be reduced below 

-58 
2 to achieve the e asked for. A 

single precision routine cannot do 

this. Reduce e and try again. 
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Information Available Upon Exit From This Routine: 



Locations 
a+i 

a+n+i' 

a+2n+i 
a+3n+i 

a+4n+i 

a+5n+i 

?8L, 80L 



Contain rounded values of variables at end of just 

completed step. 

Contain approximations to scaled derivatives ;2~ f 

(not best values) pertaining to end of just completed 

step. 

Contain carried-over truncation information. 

Contain rounded values of variables pertaining ifco 

beginning of just completed step. 

— m — 
Contain the averaged scaled derivatives 2 f 

(best slope of chord ) over interval just completed . 

Contain the 3rd approximations to 2" f , against 

which the N(a+4n+i) were tested to establish that the 

interval was adequately small. 

Contain the right address >*--m, specifying the interval 

JLast chosen for use. 



Notes On Scaling and Choice of m, iL, e: 

In general if less than full register precision is required, it is best 
to scale the y.'s so that they are small enough so that m can be taken small. This 
in turn allows Jl~ to be taken small if desired for increased speed. Truncation 
errors will not propagate more than a few places into the y. registers because 
sucft errors are taken into account. 

X should be chosen as small as possible for speed, consistent with 
printout or display intervals desired. If some 2" f . has a relatively sharp 
fluctuation in value in a narrow range of x, then one must guarantee, by^ taking 
J? large enough or otherwise, that at least ,ne point x = (integer) 2 
falls within the region of fluctuation, otherwise there is some risk of the 
program missing the fluctuation entirely. 



e should be chosen after the scale for the y. has been chosen., If the 
maximum value of |y. | is 2 } then there will be between e-d and [3A e]-d 
significant bits of the y. developed (for Ax = l)^ depending on how ill-behaved 
the equations are. As a further guide in the choice of e? we may remark that for 



extremely ill-behaved equations the interval will occasionally be forced down to 
■[3.' 

■eA 



-[3A el-m 
but not beyond 2 ' • For well-behaved equations the interval is only forced 



down to about 2 

The auxiliary must be written to provide precision consistent with the 
precision demanded of this routine ^ namely e bits or 39 bits^ whichever is the 

lesser, correct in 2 f . . 

s 1 

Method of Finding a Root of F(x s y. (x)) = 0: 

1. Integrate until F reverses sign. 

2. Back off by copying H'£a + 3n + i) back into a + i; read JC and put 

Ji_ equal to jf + 2 or 3 (3 is probably the best) 5 re-enter subroutine 
at right side of 5L» 

3. When Jc has increased to or beyond some chosen final value < [3A e] - 
1 + m and < 38., terminate loop and read root x from location a = S3. 

DESCRIPTION OF SUBROUTINES 

This program is designed to speed up the integration of ordinary 
differential equations (or of definite integrals) by choosing automatically, and 
independently for each integration step^ a near optimum value of the "interval" s 
i.e. the largest value of increment of independent variable consistent with the 
accuracy desired. The choice of interval* is limited to inverse powers of 2 in order 

to make the arithmetic 'fast in* a binary machine,, but intervals all the way from 

=2 ~^8 
2 to 2 " may be chosen. Thus if the optimum interval varies widely wer the 

complete range of integration much time is saved compared to a similar integration 

in which the minimum interval is used throughout , 

Whether or not the optimum interval varies widely 5 the user is also largely 

relieved of the necessity of estimating the correct interval for the accuracy desired, 
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The "price" for the above described automatic interval feature is 2n 
additional temporary storage locations and an average of 1 discarded integration 
step for every k useful steps. 

The routine, upon being entered finds an optimum interval and executes 
one step at this interval and returns control to the link address. Thus in order 
to perform a sequence of integration steps a master routine must be employed and 
this master routine should monitor the independent variable, not a step counter, 
to determine how far the intejkration has proceeded. 

We may integrate an# set of n simultaneous flr$b order ordinary 
differential equations (within the memory limitations of the computer) in which 
each derivative is expressed explicitly in terms of (or at least calculable from) 
the *vari4||ies themselves: 

dx , 

r- f lk' y l'>2 -■ y n-i ) 



dx 



••(MlfclWlP* * 



rf n-l 



- f H~l (X ' y l' y 2 ■— W 



dx ft 

It is necessary to include the indej»endent variable equation whether or not the 
other derivatives f to f depend explicitly upon x because the interval will 
be unpredictable and in general variable. This rules out a counter for x} however 
for speed and accuracy the x equation is handled as an integer operation. 

The computation of the functions f. to f , , given tfefcir arguments, is 
performed b^ an auxiliary subroutine which must be supplied by the user since it 
depends on his particular problem. This auxiliary reads the variables from locations 

a+ij | a 0, 1, ... n-1 and writes the scaled derivatives 2*" f. Into locations 

»m 

a + n |i, i = 1, 2, ... n-1. Note that the auxiliary need not; Supply 2 at 

location a + n. The auxiliary is called in k times in the course of any one 
integration step, whether this step is used or discarded because of failing the 
interval test. 
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Since in a sequence of integration steps the "net" change in y t may be 
the sum of many numbers (one for each step), measures Tare I'-^^^o-^^^aM^e^te 
accumulation of truncation errors in 1to j^o^s- Locations a + 2n +1, i > 0, 
are used for this purpose and hold essentially l/2 plus the least significant 
part of the 'current net increment of y i - These locations require no attention on 
the part of the user. Locations a + 1 meanwhile hold the current rounded y^* 

The standard Runge-Kutta Jrth order scheme* is used- In this scheme four 
approximations z' m f[ , 2" m ff , 2" m f* n , 2> f f , each^depending on the earlier 
ones, are successively computed and the best mean value 2~ m I ± for the interval is 
then computed as a weighted average of these four with weights l/6, l/3, l/3, l/6> 
Now of the four successive approximations, the third 2" m t ± is nearest the average 
2" m f for small enough interval and differs from it by terms of order (interval) * 
Therefore this 3rd approximation when acquired is held and when 2~ m t^ is acquired 
a test is made: 

2' U |f i .fJ II |<2- [5 ^ e]+i - m , alli>l. 

where 2"^ is the interval being tried . If this test fails the integration is 
discarded, the interval is halved and another integration tried. A failure of 
the test when i - 38 leads to an FF stop, since x can not be integrated properly 
with any interval smaller than 2~™ because the mid-point of the interval must 
be representable in the machine. If the test succeeds the integration step is 

accepted. 

The preceding paragraph describes the rules by which the interval is 
decreased when necessary. The routine also automatically increases the interval 
when possible, according to the following rules: if an acceptable integration 
step has been performed with interval 2"^ , if £<# Q , If ^ e ^ th Mt of x 
is and if a delaying binary counter is right, the interval is doubled and another 
integration is tried, subject to discard if the test fails. The reason fo l£ 
monitoring the j£ th bit of x is to insure that the points x = (integer) 2 
will not be passed over. The reason for the delaying binary counter is to prevent 
excessive numbers of waste integrations in attempting to increase the interval 



♦ 



See Milne: Numerical Solution of Differential Equations, p. 72 for a description 
of the ^th order Runge-Kutta method used. 
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too soon after it has had to he decreased. The counter plus the X-th bit of x 
cause about four integration steps to have gone by since the last interval 
decrease, before an interval increase is tried. _6 

When the routine is entered the first time it tries the invertal 2 
and uses it or halves it and tt'i es again in accordance with the outcome of the test. 
When the routine is entered by the normal entry after it has done some integrating, 
it remembers the interval last successfully used and tries that same interval first 
(or double that interval if all the conditions for attempting to double are satisfied) 



DAT E April 8, 1958 RT: IO/7/60 
PROGRAMMED B Y A. T. Nordsieek 
APPROVED BY D. E. Muller 



nj 



LOCATION 



ORDER 



3 
k 

5 
6 

7 



10 

11 

12 

13 

1^ 

15 
16 



"PP 1 

LO 5F 
50 T7L 
00 IF 

kO F 
00 IF 
iA F 
lA 3F 
L^ 63L 
42 63L 
LO 63L 
kO IF 
L3 IF 
32 TL 
FF 68F 
L5 F 
LO 7TL 
32 9L 
FF 69F 
LI 5F 
36 16L 
L5 78L 
LO 5F 
32 12L 
FF 70F 
LI 6f 
32 16L 
F5 78L 
LO 6f 
LO 5F 
36 17L 
FF 71F 
FF 72F 
FF 73F 



NOTES 



PAGE 1 



F 6 



Interlude 



Test a + 611 



Test n 



Test m > 1 



Test m < 37 



_ Test p n > m + 1 



Test t Q < 38 



LOCATION 


J : _ 

ORDER 


1 NOTES PAGE 2 


F 6 


17 


LI IF 
10 2F 


m 


Interliide 




18 


L4 7F 

fo 6f 




Test e 




19 


56 20L 
FF 7^F 








20 


09 IF 




■waste 




21 


10 S5 
40 IF 
Fl IF 




2" m -2" 59 to location (a + n) 




22 


4o S4 
F5 78L 








23 


LO 5F 
^2 Y5L- 








24 


F5 5F 
L4 f 








25 


42 71L 
00 20F 








26 


^6 71L 
LO 79L 








2? 


46 65L 
L5 4f. 








28 


L4 F 
42 56L' 








29 


42 58L 
00 20F 




Plant addresses depending ona 3 


n s m 


50 


46 64l 
46 74l 








51 


L4 79L 
46 70L. 








52 


L5 5F" 
L4 f 








55 


L4 F 
k-2 59L 









■M 



LOCATION 


ORDER 


t 


. ,11 , = ! *-l 

NOTES PAGE 3 F 6 


3h 


F5 59L 
42 60L 




Interlude 


35 


42 66L 
00 20F 






36 


46 60L 
46 66l 






37 


46 68L 
LO 79L 






38 


46 65L 
46 69L 






39 


L5 ^F 
L4 f . 






40 


L4 f 
42 57L 






41 


42 6lL 

42 67L 






42 


00 20F 
46 6il 






\y 


46 67L 
F5 67L- 






44 


42 62L- 








42 68L 






45- 


L5 55L 
4o 107L 






46' 


F5 k5L - 








L4 79L 




Copy constants into position 


vr 


4o 45L 
LO 5^L ' 






48 


36 45L , 








LI 7F ' 






h9 


10 2F 
L4 7F ■ 






50 


LO 6f 




Construct and place test quantity 




42 51L 


, ,-,, .... , 


, , ..,. ,,. — 



LOCATION 


ORDER 


NOTES PAGE 4 F6 


51 


09 IF 




Interlude 


- 


10 F 


by 50' 




52 


ko lO^L 








kO 105L 




Reset delaying counter 


53 


19 37F' 
26 80L 






5^ 


75 76L 
kO 128L 






55 


50 100L 
10 s6 






56 


N2 2L 
^9 F 






57 


N2 80L 
kl F 






58 ~ 


L5 S3 
40 F 






59 


75 Sh- 
ho F 






60 


IA F, 
40 F 






61 


74 F 
kO F 






62 


L5 1S4 
kO F, 






63 


75 F 
40 F 






6k 


lA F 

< 

40 S3 




Blank constants 


65 


74 F 
40 sV 






66 


7J F 
IA F 






67 


LJ F ' 








lA F 




.J .i , . — ■ , ... — ■ 



LOCATION 


ORDER ] t NOTES PAGE 5 


F6 


68 


I 

40 F 

LO F 


1 


Interlude 




69 


75 F. 
40 s4 








70 


L4 F 
40 IS 3 








71 


L4 f 
40 F 








72 


F5 IS 3 
40 1S3 








73 


65 S^ 
40 S4 








74 


L5 F 
40 S3 








75 


50 100L 


J 






10 F 








76 


00 F 
00 L 




Relativizer 




77 


00 F 

00 4f 








78 


00 F 
00 37F. 








79 


00 IF 
00 F 








80 


4o 128L 
50 76L 




2~ 5 to location 128 




81 


L5 71L 
42 82L 








82 


42 2L 




Reset locations a + 2n + i to l/2 




49 F 


"by 80' 






83 


F5 82L 
40 82L 








84 


LO 56L 
32 82L 




-. 




1 85 

1 


26 999F 
26 L 




Ind of Interlude 




i 


26 IN 









LOCATION 


ORDER 




NOTES 


PAGE 6 





K5 F 
42 92L 


4 


Normal ientry 




1 


26 12L 
F5 4f 




initial entry 






^ 




2 


42 2L 










49 F 


by 2,3' 


reset locations 




3 


F5 2L 
40 2L 


* 


a + 2n + i 
to 1/2 




4 


LO 108L 








5 


32 2L 
L5 6F 




interval modifying entry 










42 T8L 








6 


42 107L 
K5 33L 








7 


42 92L 
LI 7F 








8 


10 2F 
L4 7F 








9 


lo 6f 

42 10L 




Construct test quantity 




10 


09 IF 

10 F 


by 9' 






n 


40 104L 
40 105L 




reset delaying counter 




12 


L5 107L 
lo 78L 




Test 1 = l Q 




13 


32 21L 
L5 78L 








14 


l4 5f 
42 15L 




Test & bit of x 




15 


Fl S3 

OOF 


by 14' 






16 


32 21L 










L5 105L 




I Test counter 





F 6 



LOCATION 


ORDER 


NOTES 


PAGE 7 


XT 


32 18L 
kl 105L 


- 


J 

Advance counter 




18 


22 21L 
L5 7$L 








19 


LO 100L 
k2 78L 




double interval 




20 


L5 lO^L 
10 IF 








21 


40 104L 
L5 HOL 








22 


40 i4-9L 
F5 21L 




Copy N(a+i) 




23 


te 51L 
F5 23L 




to a + 3n + i 




2 V 


22 hQL 
L5 112L 




t 




25 


k2 26L 










F5 78L 




Shift address for passes 2 and 


3. 


26 


1+2 80L 
kl F 


by 25, 27 1 - 






87 


F5 2bL 
kQ 26L 




Clear locations 
a + kn + i 




£8 


LO 109L 
32 2DL 


~ 






29 


L5 101L 
J+o 106L 




Reset pass counter j 




30 


k2 39L 
50 30L 




th 
shift address for j pass 




31 


26 S6 




_. Call in auxiliary 






Ff; lOfaL 








32 


LO J.02L 
kQ F 




Test for j - 3- 


i 


1 J J 


L3 F 




i.e. 3rd pass 


i 


1 ; 


36 45L 


iii i ' 




} 



F6 



LOCATION 


ORDER 


NOTES PAGE 


8 


34 


F5 4F 
42 38L 










35 


L5 112L 












40 40L 


• 




Accumulate 3/4 2~ m f in 




56 


F5 24L 
42 43L 






locations a + 4n + i 




37 


00 IF 
F5 48L 










38 


42 44l 










39 


L5F 
L4 128L 

;of 






~J Add 2" for average round-off 




4o 


L4F 
40 F 






Closed subroutine for 




4l 


F5 38L 
40 38L 






X 2"^ + Y -» Z 




42 


L5 40L 

L4 iool 


. 








43 


40 40L 
LOF 










44 


32 38L 










45 


*22F 
L5 78L 




— 1. 


Shift address for pass 4 






42 80L 




1 
-J 




46 


L5 H4L 
40 4% 






Copy 2" f . to locations 




47 


L5 82L 
42 51L 






a + 5n + i for test 




48 


L5 6l 
42 52L 










49 


L5F 
40F 










50 


L5 49L 
L4 IOOL 






Closed subroutine for copying 





F 6 



LOCATION 



ORDER 



51 


40 49L 




LO P 


52 


56 49L 




22 p 


55 


FF 75F 




F5 106L 


5^ 


4o io6l 




F0 102L 


55 


32 62L 




L5 4f 


56 


42 38L 




L5 80L 


57 


42 39L 




L5 h6l 


58 


40 40L 




F5 57L 


59 


42 43L 




F5 59L 


60 


26 38L 




L5 106l 


6l 


10 IF 




26 30L 


62 


FF 76F 




L5 118L 


63 


40 65L 




L5 120L 


64 


4o 66l 




50 103L 


65 


7J F 




L4 F 


66 


40 P 




LO F 


67 


40 F 




L7 F 



NOTES 



PACE 9 



T 6 



test quantity has overflowed 

Advance and test for 4 passes completed 



Next intermediate values of x. y. 

j 

to locations a + i 



J 



Shift address for passes 2, 3, 4 



1 by 63, 71 f 



by 64, 70 



Interval too small 



2~ m f from 
3/4 2" m f and 

m TTT 

compare with g~ f 



LOCATION 


ORDER 


NOTES PAGE 10 


68 


IA 104L 
36 93L 








69 


L5 66L 
L4 100L 








TO 


4o 66L 
L5 65L 








71 


L4 100L 
40 65L 


- 






72 


LO 119L 








73 


32 64L 
L5 112L 
42 77L 




— 




lh 


L5 122L 
40 79L 








75 


L5 123L 
40 8lL 








76 


42 83L 
L5 124L 






Double precision addition 


77 


40 84l 
L5 F 


by 73', 85 » 




of increment 


78 


50 100L 

10 S6 


by 5 , ,19 I ,96* 






79 


L4 F 
40 F 


"\ by 7^, 87 






! 8o 


32 80L 
S5 F 


waste 

by 26, 45 T 






81 


L4 F 

40 F 


) by 75S 88' 






* 


36 85L 
36 115L 








03 


L4 10] L 






y 


| 

84 


40 f 
F5 F 


oy 16, 89 

i 




> Carry into a 




40 F 


J by 77, 90' 





J 



T6 



LOC ATION 

86 

87 
88 

89 
90 
91 
92 

?3 

94 

95 

96 

97 
98 



100 



101 



I* I i 



ORDER 



F5 77L 
40 77L 

L5 79L 
L 1 * JOJL 
4o 79L 
L5 8lL 
L4 100L 
ko 81L 
42 83L 
L5 84l 
L4 100L 
40 84L 
LO 125L 
52 77L 
00 IF 
22 F 
L5 104L 
50 101L 
00 IF 
40 104L 
36 53L 
40 105L 
F5 78L 
40 78L 
FO 12 7L 
36 62L 
L5 126L 
40 4-9L 
F5 63L 
26 2% 
00 IF 
00 IF 
80 F 
00 y )V 



NOTES 



PAGE 11 



P6 



Waste 
Exit link 



Double and test the test quantity 



Reset counter 



Halve and test interval 



Copy back x, y for new 
t r ia 1 int egra t i on 



Address increment 



Starting constant for pass 1 



LOCATION 



102 



103 



104 
105 
106 
107 

108 

109 
110 
111 
112 
113 
111*. 
115 



ORDER 

80 F 

00 6F 

00 F 

00 1832 5193 7963F 
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P6 



Ql 129K 



-2 



-[5Ae]+i2-iB 



counter 

-1+0+2 )2 -59 

50 100L 
10 A -m 
N2 2L 
k9 a+3n 
N2 80L 
iH a+5n 
L5 a 
kO a+3n 
75 a+n 
kO &+kn 
lA a+tai+1 
kO a+4n+l 
jk a+5n 
kO a+5n 
L5 a+n+1 
kO a+5n+l 
75 a+2n 
kO a+6n 



Test constant for pass 3 and end 



1/3 



Causes SADOI to put following material 
on problem tape immediately after this 
program* 

Test quantity- 
Used to delay interval doubling 
j=l, . .. k for passes 1 ... k 



starting, end and test constants 
formed by interlude 



LOCATION 


ORDER 
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F6 


116 


L4 a+3n 
40 a 








117 


74 a+4n 
40 a+n 








118 


7J a+4n+l 
L4 a+4n+l 








119 


LJ a+5n 
L4 a+5n 








120 


40 a+4n+l 
LO a+5n+l 




Starting, end and test constants 




121 


75 a+4n 
40 a+n 




formed by interlude. 




122 


L4 a+3n+l 
40 a+1 








123 


14 a+2n+i 
40 a+2n+l 








124 


F5 a+1 
40 a+1 








j.25 


65 a+n 
40 a+n 








126 


L5 a+3n 
40 a 








127 


50 100L 








128 


10 38-m 
2 -58 




. 








. —4- . . , . 





